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Abstract: The atom-to-molecule conversion by the technique of optical Feshbach resonance in a magnetic 
lattice is studied in the mean-field approximation. For the case of shallow lattice, we give the dependence of the 
atom-to-molecule conversion efficiency on the tunnelling strength and the atomic interaction by taking a double- 
. well as an example. We find that one can obtain a high atom-to-molecule conversion by tuning the tunnelling 

■ and interaction strengths of the system. For the case of deep lattice, we show that the existence of lattice can 

| improve the atom-to-molecule conversion for certain initial states. 
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I. INTRODUCTION 



Since Bose-Einstein condensations (BECs) in dilute atomic gases were realized in 1995, the study of cold atoms becomes 
a remarkable research area which has been extended from atomic to molecular systems in recent years. Molecular BECs are 
versatile not only for cold atomic physics but also for other research areas (H-d], because they include more degrees of freedom 
than atomic systems. To realize molecular BECs, one usually converts ultracold atoms into molecules through resonant photoas- 
sociation (optical Feshbach resonance) or magnetoassociation (magnetic Feshbach resonance) Jj-Ql]] rather than cool molecules 
<4— > . directly. Additionally, besides the technique of Feshbach resonance, the stimulated Raman adiabatic passage technique was also 
used in the conversion of Bose-Fermi mixture into molecules lfl2l[T3ll . In the above works, the atom-to-molecule conversion 
system were all confined in a single well. Note that in recent experiments iTT^l — TTvIl . the atom-to-molecule conversion system 
confined in an optical lattice was also studied, where the atom-to-molecule conversion efficiency and the lifetime of molecules 
can be improved due to the suppression of the inelastic collisions. Similar to optical lattice, magnetic lattice can also be expected 



to improve the atom-to-molecule conversion since ultracold atoms have been successfully transferred to the magnetic lattice po- 
tential experimentally in 2008 [18]. Compared with optical lattice, magnetic lattice has several distinct advantages, such as high 
stability with low technical noise and low heating rates, large and controllable barrier heights and so on. It is thus meaningful to 
study the property of atom-to-molecule conversion in a magnetic lattice. 

In this paper, we consider an atom-to-molecule conversion system in a magnetic lattice. In the mean field approximation, 
we study the influence of magnetic lattice on the atom-to-molecule conversion efficiency and show how to improve the atom- 
| to-molecule conversion in a magnetic lattice. The rest of this paper is organized as follows. In Sect. [II] we present the model 
I/-) ■ and the dynamical equations. In Sect. |IIT] we consider the case of shallow lattice, i.e., the atoms can tunnel between the nearest 
neighbouring sites. Taking a double-well as an example, we study the time evolution of molecular density and the effect of 
parameters of the system on the atom-to-molecule conversion. We also confirm our numerical results with the help of fixed 
points and energy contours of the system. In Sect. [IV] we consider deep lattice, i.e., the atoms can not tunnel between the lattice 
sites. We show that the existence of the lattice can improve the atom-to-molecule conversion for certain initial states. In the last 
section, we give a brief summary and discussion. 



H. MODEL AND GENERAL FORMULATION 

We consider the atom-to-molecule conversion via optical Feshbach resonance in the magnetic lattice, where the atoms are 
subject to the lattice but the molecules are not. The Hamiltonian describing such a system can be written as, 

j. -t. U a v~i U b g v~i ~ ~ + + ~ -i- ~ 

Hi = -J 2J< a i a i + a j a >) + "jj - "<"'("<"' ~ ^ + —h{n b -\) + -2_ u (b'a i a i + ba\a\) + 5b f b, (1) 

where the operator al^ and creates and annihilates a bosonic atom in the i(y')th site separately and (z, j) denotes the two 
nearest neighbouring lattice sites, while ¥ and b creates and annihilates a bosonic molecule. They obey the commutation relation 
[a,-, fit] = 6ij and [b, b~' f ] — 1. Here the parameter J denotes the tunneling strength of atoms between the nearest neighbouring 
sites, which can be tuned by changing the distance between the neighbouring lattice sites or the height of the potential barrier 
separating the neighbouring lattice sites. U a and Ui, is the interaction strength between the atoms and molecules, respectively. 
As we know, U a can be written as U a - Intp-alu, where a and fi are scattering length and reduced mass, respectively. The 
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scattering length can be tuned by the optical Feshbach resonance technique. 6 is the energy detuning between the atomic and 
molecular states, and g refers to the coupling strength between atoms and molecules. The molecules in spin singlet state (S = 0) 
are not subject to the magnetic lattice, so there is no subscript for b in the Hamiltonian (HJ. 

In order to study the property of atom-to-molecule conversion, we need to give the dynamical equations of the system. The 
classical field description is an excellent approximation if the quantum fluctuation is small. As we know, the magnitude of 
quantum fluctuation around the condensate state scales down as 1 / VA 7 in zero temperature with N the particle number. There 
are usually 10 4 - 10 7 particles in dilute BEC experiments, so we adopt the mean-field description 11911 . In the mean field 
approximation, by replacing operators with their expectation values, i.e., a^-j — > (S;(n) = auj> and b — » (b) — b, one can easily 
give the dynamical equations for 2,-yj and b with the help of Heisenberg equations of motion for and b. These equations 
obey the conservation law 2/ aj&i + 2Pb = N. To simplify the calculation, one usually assumes a^j) = a^/ VA 7 and b — bj VA 7 
with 2, |fl/| 2 + 2|b| 2 = 1. Then the dynamical equations for a%j) and b can be rewritten as, 

idi — -J(di+\ + fl/_i) + U a \ai\ 2 cii + gba*, 

A = U b \b\ 2 b + 6b + ^Yj a l ( 2 ) 

i 

where U a = U a N, Ub = UbN, g = g VA 7 and natural units has been used. Note that the interaction strength between molecules 
is much smaller than that between atoms in most experiments, so we ignore the interaction between molecules in the following 
calculation, i.e., U/, = 0. Meanwhile the energy detuning is chosen as 6 — . We study the atom-to-molecule conversion 
efficiency, which is defined as twice the density of the largest molecules (2\b\f mx ) in the time evolution, for the cases of shallow 
lattice and deep lattice in the following two sections separately. 

III. SHALLOW LATTICE 

In this section, we consider that the atoms tunnel between the nearest neighbouring sites, i.e., J + 0. Although we can give 
the dynamical properties of the system for different number of lattice size L by solving numerically, here we mainly take 
double-well as a typical example. For the case of double-well, Eqs. is simplified to, 

ia\ — —Jci2 + U a \ai\ 2 ai + gba\, 

ici2 = -Jci\ + Ua\a2\ 2 a2 + gba* 2 , (3) 
ib = g{a\ + a%)l2, 

where 6 — and Ub = have been taken. Tunnelling strength /, atomic interaction strength U a and the atom-molecule coupling 
strength g have the same dimension. We choose g as unity and then all quantities are renormalized to be dimensionless. 

A. Evolution and Conversion Efficiency 

Now we are in the position to study the atom-to-molecule conversion of the BECs in double wells with the help of (0. Here 
we consider the symmetrical initial state a\ 2 = Vl/2 and b — 0. We plot the time dependence of molecular density for different 
parameters in Fig.Q](a), (b) and the period of Rabi oscillation [20] versus U a in Fig.[T|(c), (d). Fig. [U(a) shows that all atoms 
can be converted into molecules and the system will always stay in the pure molecular state, i.e., the oscillation period for this 
condition is infinite which corresponds to the peak in Fig.Q](c). Whereas, from Fig.Q](b), not all of atoms can be converted 
into molecules for some parameters and the system oscillates periodically are found. The oscillation period for this condition is 
finite corresponding with the peak in Fig.[T](d). The above results enlighten us on choosing the appropriate parameters for high 
atom-to-molecule conversion efficiency. So it is necessary to find the suitable relation between tunnelling strength J and atomic 
interaction strength U a - 

We plot the dependence of atom-to-molecule conversion efficiency on the atomic interaction strength in Fig.|2](a) and (b), and 
on the tunnelling strength in Fig.|2](c) and (d). Fig. [2] (a) and (b) show that for 7^0, the increase of interaction strength can 
improve the atom-to-molecule conversion when U a is smaller than the critical value U ac , while suppress the atom-to-molecule 
conversion when U a > U ac , which is different from the case of J — plotted in Fig. [5] (b). Additionally, for a given value 
of J, when U a = U ac , one can get the highest atom-to-molecule conversion efficiency which is dependent on the value of 
J. Comparing Fig. [T](c) and (d) with Fig. [2] (a) and (b), we find that both the periods of Rabi oscillation versus U a and the 
dependence of the maximum atom-to-molecule conversion efficiency on U a exhibit a similar behaviour and they have the same 
parameters of the peaks, e.g., J = 0.5g, U a - 2g for Fig.Q](c) and Fig.|2](a) and J — g, U a — 3.853g for Fig.Q](d) and Fig. [2] 
(b). When J = Q.5g and U a = 2g, all of atoms can be converted into molecules and stay in the pure molecular state finally with 



3 




U a (units of g) W U a (units of g) (d) 

FIG. 1: (Colour online) Time dependence of the molecular density with symmetrical initial conditions in (a) and (b). The parameters are 
J = 0.5g and U a = 2g (a), J = g, U a = 3.853g (b). The period of Rabi oscillation T versus U a in (c) and (d). The parameters are J = 0.5g (c) 
and J = g (d). 




J (units of g) ( c ) J (units of g) ( d ) 



FIG. 2: (Colour online) Dependence of the atom-to-molecule conversion efficiency on the atomic interaction strength (a) , (b) and the tunnelling 
strength (c), (d) for double-well. The parameters are J = 0.5g (a), J = g (b), U a = 2g(c), and U a = 4g(d). The vertical pink dashed lines are 
guides to the eye. 
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infinite period. While only partial atoms can be converted into molecules and the period is finite with parameters 7 — g and 
U a = 3.853g. Fig.|2](c) and (d) tell us that for the case of U a + 0, the dependence of atom-to-molecule conversion efficiency 
on the tunnelling strength does not change monotonically but there is a critical value J c where the atom-to-molecule conversion 
efficiency is the highest. 




FIG. 3: (Colour online) Projection of the atom-to-molecule conversion efficiency in the J-U a plane for the cases of double- well (a) and triple- 
well (b). The red solid line and red dotted line distinguish the complete and uncomplete conversion respectively. The purple filled points are 
the general critical point. The black dashed lines are U a = 4 J (a) and U a = \2J (b) which are guides to the eye. 

In order to give further overall dependence of atom-to-molecule conversion on the strengths of tunnelling and interaction of 
the system, we plot the projection of the atom-to-molecule conversion efficiency in the J-U a plane in Fig. 0(a). The purple filled 
point {J'c'Uac) = igl V2,2 V2g) in Fig.[3](a) is an important critical point called "general critical point". When 7 < $ c , all of 
atoms can be converted into molecules if the parameters are selected as the red solid line, i.e., U a = 47. On the other hand, for 
the case of 7 > J~ c , although not all of atoms can be converted into molecules no matter what parameters are taken, one can also 
get higher atom-to-molecule conversion by taking the parameters on the red dotted line which deviates from the line U a = 47 
(i.e., the black dashed line). Fig.[3]shows both the atomic interaction and the tunnelling strength can affect the atom-to-molecule 
conversion efficiency and one can get high conversion by tuning the parameters of the systems properly. In order to convert all 
atoms into molecules, one should choose the parameters on the red solid line, i.e., U a = 47 and 7 < JT C - 

Additionally, we extend our study to triple-well case and give the projection of the atom-to-molecule conversion efficiency in 
the J-U a plane in Fig.[3](b). There is also a general critical point (J'cUac) = (g v2/4, 3 V2g) and the relation of the red solid 
line is U a - 127 as 7 < JT C . The discussion of the dependence of atom-to-molecule conversion efficiency on U a , J and the 
oscillated period are similar to the double-well case. 
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B. Fixed Points and Energy Analysis 

In order to understand the dependences of the conversion efficiency on the parameters J and U a , we study the property of the 
fixed points. By expressing a,- = -\fp a le e ' ,i and b - -\fphe 6 '', we can write Eq. @ as 

z = -27 ^PaiPai sin(20 fl ) - g yfp~b[p a i sin(0 - 2<p a ) - p al sin(0 + 2cj> a )], 

COS(20 fl ) + | -y^[cOS(0 - 24> a ) - COS(0 + 2(f> a )] + ^Z, 



yjPa\Pa2 2 2 

2/(1 -2p b ) 



Pb = -y/pb[p a ism((f)-2(f> a )+p a 2sm((f) + 2^ a )], (4) 



cos(20 fl ) - g -^[cos(0 - 2<p a ) + cos(</> + 20 fl )] - t/ fl (l - 2pi) 



VPalPa2 

Q 

■_[p a i cos(0 - 20 a ) + p a2 cos((/> + 24> a )] 



(5) 

where z = p a \ - p«2, <Pa = (0 a 2 - S a i)/2, 4> = 6 a] + 6 a2 - b and p al = (1 - 2p b + z)/2, p fl2 = (1 - 2p fe - z)/2. z and a , p 6 
and (p are mutually canonical conjugations, respectively. They satisfy the Hamiltonian canonical equation, i.e., z = -dH c i/d(f> a , 
<p a = dH c ildz, Pb = -dH c i/d</> and (f> = dH c ijdp b . The corresponding classical Hamiltonian of the system is given by 

H cl = -JiJ(l- 2p b ) 2 - z 2 cos(20 fl ) + -^[(1 - 2 Pi ) 2 + z 2 ] + | ylpb[p a i cos(<f> - 2<j) a ) + p a2 cos(0 + 20 fl )]. (6) 

We should note that the time evolution of system can be determined by the classical Hamiltonian due to the energy conservation 
law. 

In order to get the fixed point solutions, we set z = 0, <p a — 0, p b = and = in Eq. For the initial states z — 0, <p a = 0, 
the atomic distribution and phase in the double wells are the same at any time, i.e., z(t) = and (f>Jf) — 0. Then the conditions 
for the fixed points can be written as, 

g^fPbiX ~ 2p b )sm(f> = 0, 
2J -U a (l-2p b ) + S \ J;"' cos4> = 0. (7) 

Here we do not write the fixed point solutions in order to save space, although we can get such solutions through solving Eq. (0. 
Note that when < J < gj y2, there is a fixed point p b = 1/2 with <f> being not well defined. Since the fixed point p b = 1/2 
exists for the case of < J < gj V2, it is possible for the system to reach it, i.e., all of atoms can be converted into molecules 
by tuning the parameters of the system properly. According to Eq. ©, we know that the energy of the system is -J + U a /4 for 
the initial state we considered, and zero for the fixed point p b — 1 /2. So if the system can reach the fixed point p b — 1 /2, the 
parameters must satisfy -J + U„/4 = due to the existence of energy conservation law. It explains well why the conversation 
efficiency corresponding to the parameters taken on the red line in Fig. [3] is 1. Whereas, for the case of J > gj y2, the fixed 
point p b - 1/2 disappears, and then not all of atoms can be converted into molecules no matter what atomic interaction strength 
is taken, which is confirmed by Fig.|2](b). 

In Fig. [4] we plot the energy contours in the phase space of p b and <f> for the two cases of J < gj V2 and J > gj y2, respectively. 
From Fig.Q](a) and (b), we can see that, for a fixed value of J, one can get highest atom-to-molecule conversion efficiency when 
U a = U ac - So in Fig. [4] we take the atomic interaction strength as U ac which is determined by J. Note that we can obtain 
Vac — 4 J analytically for J < gl y2, while the value of U ac needs to be determined by numerical method for J > gj y2. For the 
initial states we considered, the system evolves along the pink dashed curve. Fig.|4](al) and (bl) further confirm that the system 
can reach the state p b = 1 /2 for J < g/ V2 but can not for J > g/ V2 clearly. 

We also study the property of the fixed points to confirm the dependence of the conversion efficiency on the parameters J 
and U a for the triple-well from Fig. [3(b) and they are in good agreement. The fixed point p b = 1/2 exists for the case of 
< J < g V2/4 when the atoms can be converted into molecules completely with corresponding critical interaction strength, 
while the fixed point p b = 1/2 disappears for the case of J > g V2/4 when only part of the atoms can be combined to molecules. 



IV. DEEP LATTICE 



In this section, we consider the case of deep lattice, where the height of potential barrier separating the neighbouring sites 
is so high that we can ignore the tunnelling of atoms between neighbouring site, i.e., J — 0. For a deep lattice without atom- 
to-molecule conversion, the ground state has a fixed number of particles per lattice site, and the relative phase between lattice 
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FIG. 4: (Colour online) Energy contour in the phase space of t/> and p b . The parameters are / = 0.707 and U a = 2.828 (a), 7 = 0.708 and 
U a = 2.832 (b). The enlargements of the small blue boxes in (a) and (b) are (al) and (bl), respectively. The pink dashed curves correspond to 
the energy with the initial states a\ = ai = 1/ Viand/? = 0. 



sites is smeared out. However, the existence of atom-to-molecule conversion, where the molecule are not subject to the lattice, 
makes the phases between lattice sites related and well defined. In order to study the atom-to-molecule conversion efficiency, 
we assume there is no molecule in the system at the initial time and the atoms are equally populated in the lattice, i.e., the initial 
state is b — and a, = y/N/L with L being the number of lattice. Once the initial state is fixed, we can get the time evolution 
of the system through solving Eq. (fJJ numerically and the dependence of the atom-to-molecule conversion efficiency on the 
atomic interaction strength U a if the value of L is not very large. We summarize our results in Fig. 0(a) and (b) for the case of 
1 = 1,2,3,4. 

In Fig. 0(a), we plot the time dependence of molecular density for different L with certain atomic interaction strength (U a = 
2g), and plot the corresponding dependence of the atom-to-molecule conversion efficiency on the atomic interaction strength U a 
in Fig. |5](b). The molecular density oscillating periodically with time for any value of L is displayed in Fig. 0(a). The larger the 
number of lattice site, the longer the oscillation period. From Fig. 0(b), we find that the atom-to-molecule conversion efficiency 
becomes smaller when the atomic interaction strength increases. For the same atomic interaction strength, the atom-to-molecule 
conversion efficiency becomes higher with the increase of the number of lattice site, which shows that one can improve the 
atom-to-molecule conversion by magnetic lattice even if the atomic interaction strength is fixed. 

In order to confirm the above numerical results, we study the system with analytical method. For the symmetrical initial state 
a, ■ — 1/ VZ and b — 0, the lattice sites are equivalent, i.e., the values of a, at any time does not change with different ;. Then we 
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FIG. 5: (Colour online) Time dependence of molecular density for different values of L with symmetrical initial state (a) and asymmetrical 
initial state (c), dependence of the atom-to-molecule conversion efficiency on the atomic interaction strength U a for different L with symmetrical 
initial state (b) and asymmetrical initial state (d). The parameters are U a = 2g and / = (a) and (c), J = (b) and (d). The asymmetrical initial 
conditions considered are |ai(0)| 2 = 1/4, |a 2 (0)| 2 = 3/4, \b(0)\ 2 = for L = 2 and |a,(0)| 2 = 1/4, MO)! 2 = 1/4, |a 3 (0)| 2 = 1/2, \b(0)\ 2 = for 
L = 3. 



can introduce A - a, VZ. Substituting this definition into Eq. (0, we can give the dynamical equation for A and b, 

iA = l j-\A\ 2 A + gbA\ 

ib = -A 2 , (8) 
2 

with the particle conservation law \A\ 2 + 2\b\ 2 = 1. Eq. ([8]) is similar to the dynamical equation for the atom-to-molecule 
conversion system in single well except for the rescale of the atomic interaction strength. From Eq. ([8), we can see that equally 
distributing the same atomic BEC into magnetic lattice can reduce the effective interaction strength between atoms for the case 
of J = through comparing with the case of single well. Additionally, we know that the existence of the atomic interaction can 
suppress the atom-to-molecule conversion, which can be confirmed by Fig.|5](b). So for the same atomic interaction strength, 
the existence of deep magnetic lattice can improve the atom-to-molecule conversion. Because the effective atomic interaction 
strength is t/ fl /L£the larger the number of lattice site is, the higher the atom-to-molecule conversion efficiency will be. We only 
consider the symmetrical initial state above. For the asymmetrical initial state, i.e., the initial atomic distributions |a,(0)| 2 are not 
the same, the existence of deep magnetic lattice can also improve the atom-to-molecule conversion if the initial phases of a, are 
the same. Time dependence of molecule density for different values of L and dependence of the atom-to-molecule conversion 
efficiency on the atomic interaction strength with asymmetrical initial state are shown with Fig.|5](c) and (d). The asymmetrical 
initial states |fli(0)| 2 = 1/4, |a 2 (0)| 2 = 3/4, \b(0)\ 2 = for double-well and |fli(0)| 2 = 1/4, |a 2 (0)| 2 = 1/4, |<a 3 (0)| 2 = 1/2, 
\b(0)\ 2 = for triple-well are considered. However, for the case of initial a, with different phases, lattice can not always improve 
the atom-to-molecule conversion. We should note that the phases of the condensates between lattice sites are well defined in 
the Mott insulator state for the existence of atom-to-molecule conversion. The phases of the condensates between lattice sites 
remain unchanged if they are the same initially. However, the phases of the atomic condensates between lattice sites will evolve 
over time if they are different initially. 
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V. SUMMARY AND DISCUSSION 

In this paper, we have studied the atom-to-molecule conversion in a magnetic lattice. For shallow lattice, where the atomic 
tunnelling strength J is unneglected, we studied the effect of tunnelling and interaction strengths of the system on the atom- 
to-molecule conversion by taking double-well as an example. We gave the dependence of the atom-to-molecule conversion 
efficiency on the tunnelling strength of atoms and the atomic interaction for the cases of double-well and triple-well. The general 
critical points for both cases were found. We also showed that atoms could be converted into molecules completely if one 
chose appropriate parameters, i.e., the parameters on the red solid lines in Fig. [3] (a) and (b) while not completely with other 
parameters. The analyses of fixed points and the energy contour were given to confirm our results. For deep lattice, where the 
atomic tunnelling strength J is neglected, we showed that if the initial phases of BECs in different lattice sites were equal, the 
existence of lattice site improved the atom-to-molecule conversion. Considering a symmetrical initial state, we showed that the 
larger the number of lattice site is, the higher the atom-to-molecule conversion efficiency could be reached. We also confirmed 
our results by an analytical method with the help of the redefinition of a,. In a word, we gave some suggestions on how to obtain 
higher conversion efficiency with suitable relation between tunnelling and interaction strengths. 

The work is supported by NSFC Grant No. 1106741 17, No. 1 1074216 and partially by PCSIRT Grant No. IRT0754. 
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